Correlations in interacting systems with a network topology 
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We study pair correlations in cooperative systems placed on complex networks. We show that 
usually in these systems, the correlations between two interacting objects (e.g., spins), separated 
by a distance £, decay, on average, faster than l/(£ze). Here Zt is the mean number of the £- 
th nearest neighbors of a vertex in a network. This behavior, in particular, leads to a dramatic 
weakening of correlations between second and more distant neighbors on networks with fat-tailed 
degree distributions, which have a divergent number Z2 in the infinite network limit. In this case, only 
the pair correlations between the nearest neighbors are observable. We obtain the pair correlation 
function of the Ising model on a complex network and also derive our results in the framework of a 
phenomenological approach. 
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I. INTRODUCTION. 

The generic features of real-world networks (the Inter- 
net, the WWW, biological, social and economical net- 
works and many others) are a complex organization 
of their connections and the small-world phenomenon 
B II H II II H II In the networks with the small- 
world effect, a mean intervertex distance grows with a 
total number of vertices, N, slower than any power of N, 
e.g., grows as IniV. Concerning economical and social 
networks, the small- world property is often considered 
as an evidence for growing close interrelations and glob- 
alization. 

Many of the real-world networks are formed by inter- 
acting objects and demonstrate complicated dynamics. 
The dependence of the correlations between interacting 
objects in systems with a network topology on time and 
distance provides a useful information about the network 
dynamics. In the present paper we discuss general prop- 
erties of the correlations between a pair of interacting 
objects on a complex network. Recall that in interacting 
systems on lattices, with a few exceptions, pair correla- 
tions decrease exponentially with distance apart from a 
critical point, where the decrease is power-law. From a 
naive point of view, one might expect that the small- 
world property of a network would enhance pair cor- 
relations between distant objects in comparison to lat- 
tices. However, it is not the case. We demonstrate that 
correlations between ^-th nearest neighbors, on average, 
decay with I as \j(£zg) or faster. In networks, where 
the mean intervertex separation 1(N) ~ IniV, this cor- 
responds to the exponential decay of correlations with £ 
(even at the critical point). However, in networks with a 
divergent mean number of the second nearest neighbors 



[where 1 (iV) grows slower than In N) , we observe a dra- 
matic weakening of pair correlations between the second 
and more distant nearest neighbors. In these networks, 
only the nearest neighbors are strongly correlated, while 
the correlations between more distant vertices are sup- 
pressed and approach zero in the infinite network limit. 

In Section^ we consider pair correlations in an inter- 
acting system defined on the top of a complex network 
in the framework of a phenomenological approach. This 
approach allows us to understand general properties of 
pair correlations irrespective of details of an interacting 
system and the nature of interactions. In Section lllll wc 
support these results by calculations of the static pair 
correlation function of the Ising model on a complex net- 
work (more precisely, the configuration model of a net- 
work a). 



II. PHENOMENOLOGICAL APPROACH 

Let a quantity Xi(t) describe a dynamic process on 
a network, where indices i label vertices, and t is time. 
Xi(t) = Xi(t) — (Xi) describes fluctuations around the 
average value (Xi). If the system under consideration is 
in an equilibrium state, then pair correlations between 
two arbitrary vertices i and j may be characterized by 
the following correlation function: 

to 

G tj (h,t 2 ) = t^ 1 / x i (t + t 1 )x j (t + t 2 )dt 



(1) 
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In the present section the brackets (...) mean the aver- 
age over the observation time to which must be much 
larger than the maximum relaxation time in the system 
under consideration. In the equilibrium state we have 
Gij(t\,t2) — Gij(t\ — £2) in the limit to — > 00. Assuming 
the Hamiltonian dynamics, one can introduce a general- 
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ized field Hi(t) conjugated with Xi(t). The function 



Xiiih-ti) = ^o 1 dx l {t + t 1 )/dH J {t + t 2 )dt 



(dXiitJ/dHjfo)) 



(2) 



is a generalized non-local susceptibility which character- 
izes the averaged response of Xi(t\) at time t\ on a field 
applied at a vertex j at time t% . There is a simple rela- 
tionship between Gij{t) and Xij{t)'- 



Gij(t) oc Xij(t), 



(3) 



where coefficient of the proportionality is unimportant 
for our purpose. 

In general case, the total susceptibility x(t) of a system 
(per vertex) is finite in the limit TV — » oo where TV is the 
number of vertices in a network: 



x(h-h) 



1 N 

-y 

TV ^ 

»=i 



dH 3 {t 2 ) 



0(1). 



H1 — H2 — ■■■ — H 



(4) 

Here, a statement A = 0(1) means that a quantity A is 
finite in the limit TV — > 00. The susceptibility x(i) may 
diverge only at a critical point. We rewrite Eq. Q as 
follows: 



x(t) 



TV ^ 



(5) 



where the first sum in the parentheses is a sum over 
the nearest neighbors of a vertex i being at the distance 
£ij = 1 . The second sum is over the second nearest neigh- 
bors, £ij — 2, and so on. One can introduce the average 
value x{t, £), that is the average value of the susceptibility 
Xij{t) a t &ij = £ over all possible network configurations. 
In the large network limit, this quantity coincides with 
the following expression calculated for a single network: 



x(t,£) 



E 



>Xij{t) 



E 



N Zl 



(6) 



where the sums are over all pairs of vertices at the in- 
tervertex distance £, ze = TV -1 ..^ ,_ t 1 is the mean 
number of ^-th nearest neighbors. Consequently, 



x(t) = Y l z eX (t,e), 



(7) 



where zq = 1 and x(t, 0) is the average local susceptibil- 
ity. Let us first assume that all x(^) > 0- I n order to 
satisfy the condition %(t) = 0(1), the terms of the sum 
in Eq. @ must decay with £ faster than This leads 
to the following restriction: 



(8) 



It is important that here is a function of the network 
size TV, so formula JSJ) shows how the nonlocal suscepti- 
bility (and spacial correlations) vary with TV. 

If the nonlocal susceptibilities Xij {t) are °f varying 
sign, formula JHJ is not applicable. In this case, for the 
root-mean-square with averaging over all pairs of vertices 
i,j with a given tij = £, we arrive at 



Xij [tj &ij 



1/2 



i) ' <o( z ; 1/2 ) 



(9) 



In accordance to Eq. ©, the £ dependence of the 
average correlation function is the same as for the average 
susceptibility, i.e., it is described by formulas JHJl or 0. 

In sparse networks, the mean number of the nearest 
neighbors (the mean degree) z\ = k = kP(k) is fi- 
nite. In general, in networks, where the mean interver- 
tex distance 1(N) ~ TV, the mean numbers of £-th nearest 
neighbors grow exponentially with £. In particular, in un- 
corrected networks (without degree-degree correlations 
between nearest neighbor vertices), Z£ = z\(z 2 j ' z\) l ~ l , 
where z 2 = £ k(k - l)P(fc) [13 • 

In this case, even at the 
critical point of a cooperative model, the nonlocal sus- 
ceptibility and pair correlations decrease exponentially 
rapidly with £. 

The effect is even more dramatic, if the degree distri- 
bution is fat-tailed, and its second moment diverges as 
TV — > 00. In scale-free networks, this takes place if the 7 
exponent of the degree distribution P{k) ~ fc~ 7 is equal 
or below 3. In this case, already the mean number z 2 
of the second nearest neighbors of a vertex diverges as 
TV — > 00. Consequently, according to Eq. (JSJ or Eq. JJjJ, 
in the limit of a large network, pair correlations between 
the second (and further) nearest neighbors vanish, and 
only pair correlations between the nearest neighbors are 
observable. 



III. CORRELATIONS IN THE ISING MODEL 

The simplicity of the Ising model and the fact that 
many microscopic models may be mapped to this model 
make the Ising model to be very attractive. Recent inves- 
tigations have revealed that the critical behavior of the 
Ising and Potts models on complex networks strongly dif- 
fers from the standard mean-field behavior on a regular 
lattice [H El El El El El E3 In the present sec- 
tion we analyze pair correlations in the equilibrium Ising 
model on an uncorrelated random complex network. 



A. The model 

We consider the ferromagnetic Ising model: 

n = - j s i s i - HiSi ' ( 10 ) 
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where Si — ±1, and Hi is a local magnetic field at a 
vertex i. The sum is over all nearest neighboring vertices. 
The static pair-correlation function 



(S t S 3 ) - (S^ {Sj) 



(11) 



is related to the non-local magnetic susceptibility Xij '■ 



(12) 



graph is described by 2L = J2i parameters xy. Here, 
L is the total number of edges of the graph. 

Using Eqs. (|13|l and l|14|l . the parameter Xy may be 
related to parameters Xji which characterize the edges 
(sub-trees) outgoing from the vertex j as follows : 



H, 



n 



(15) 



where (3 = 1/T and T is temperature, M$ = (Si). In 
this section (...) means the statistical average with the 
Hamiltonian TL. 

As a substrate, we use the standard model of an un- 
corrected random network — the configuration model @ . 
This is the maximally random graph with a given degree 
distribution or, as it is called in graph theory, a labelled 
random graph with a given degree sequence. It is im- 
portant that the uncorrelated networks have a locally 
tree-like structure. That is, in the large network limit, 
there are no loops in a finite neighborhood of a vertex. 
These networks may be considered as the random Bethe 
lattices, which, by definition, have no boundary. One 
should note that in contrast, the Cayley tree contains 
boundary vertices which are dead ends 



B. How to solve the Ising model on a tree-like graph 

Let us outline an effective approach to solution of coop- 
erative models on tree-like networks, which we use for ob- 
taining Xij- This method was implemented for the I sing 
and Potts models on regular and random [l2L Il9f 
Bethe lattices. Consider an arbitrary tree-like graph. 
Consider a spin Si on a vertex i of this graph with fcj 
adjacent spins Sj, j — 1, 2, . . . , fej. As any vertex in this 
graph, vertex i may be treated as a root of a tree. In 
turn, vertices j may be treated as roots of subtrees grow- 
ing from the vertex i. In order to characterize the subtree 
with the root spin Sj we introduce a quantity 



9ij{Si) 



eX P 

{S,}=±1 



0J S nS: 
(nm) 

+/3JS i S j +/3j2HnS, 



(13) 



Here the indices n and m run only over spins that belong 
to the sub-tree, including Sj. The statistical sum in Eq. 
(|13fl is taken only over the spins on the subtree. 
We introduce a quantity 



= gij(-l)/gij(+l). 



(14) 



In such a way, for each vertex i we can introduce ki 
parameters xy, j — 1,2,. It should be noted that 



Xji because x^j 



and Xji characterize different sub- 



trees. So, each edge is characterized by a pair: xy and 
Xji. In sum, the Ising model on an arbitrary tree- like 



where I is the index of the first nearest neighbors of the 
vertex j. The vertex i is excluded from the product over 
I. In other words, a vertex I is a second neighbor of the 
vertex i and the first one of the vertex j. The function 
y(H,x) in Eq. (|15f) depends on a cooperative model. For 
the Ising model, 



y(H,x) 



,{-J+H)P + e (J-H)ff s 



(16) 



2L independent parameters Xy should be found by solu- 
tion of the set of 2L equations l|15|) at a given temperature 
T and local magnetic fields Hi . For an arbitrary tree-like 
graph, these equations may be solved numerically, e.g. , 
by use of the population dynamic algorithm, see Ref. |l7j 
where this method has been applied to the Potts model 
on a tree-like graph. 

Observable thermodynamic quantities of the Ising 
model may be written as functions of the parameters 
x^ . For example, a magnetic moment Mi is given by 
the following equation: 



Mi 




(17) 



i=i 



Finally, for a random network, the resulting observables 
should be averaged over various configurations with ap- 
propriate statistical weights. 



C. Derivation of the pair correlation function 



Let us find a non-local susceptibility Xij f° r tw° arbi- 
trary vertices i and j at the distance hj = i from each 
other. On a tree-like graph there is the only shortest 
way connecting these vertices. It starts from i, then goes 
through vertices i\, «2> • • • > ii-i and ends at j. Using Eqs. 
C3 and (pjj) we get 



Xij 



dMj dxui dx lll2 

dxa^ 0xi 1 i 2 9xi 2 i^ 



9xi l 2 i e i dxi e _ 1 j 



dx 



dHi 



(18) 



If we know xy we can get Xij ■ 

First, for the purpose of comparison, we find a non- 
local susceptibility Xij of a regular Bethe lattice with a 
coordination number k. 

In a uniform magnetic field H\ = H 2 = .... = H 
all vertices and all edges in a regular Bethe lattice are 
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equivalent. Therefore, the parameters Xij are equal, i.e. 



and Eq. lfT3|) is reduced to 

x = y(H, x 1 *- 1 ) 



(19) 



This equation determines function of T and H. 

From Eqs. and (jT§j) we get 



1 DM 

k dx 



1 dy{H,x k ~ x y 



(fc-1) 



<9:r 



an 



(20) 



where in accordance with Eq. (|17() we have M = (e 2/3H - 
x k )/{e 2l3H +x k ). At zero magnetic field H = 0, Eq. lO 
takes a form: 



4/3X* 



(l + z fe ) 2 



2x k ~ 2 sinh(2 J/3) 



(e^ + , 



(21) 



This equation is valid at all temperatures. 

In the paramagnetic state, i.e. at temperatures T 
above the critical temperature T c = 2J/ln[fc/(fc — 2)] 
of the ferromagnetic phase transition, the self-consistent 
equation l|19|l has the only solution x = 1, and we get 



x (£)=/3G(l) = ptanh e (J0) 



(22) 



in agreement with the result obtained in Ref. [2(j for a 
Cayley tree in the framework of another method. Fur- 
thermore, using the fact that x^ = 1 at T > T c , one can 
show that in the paramagnetic phase, x(£) of the Ising 
model on an arbitrary tree-like graph is the same as that 
for an arbitrary Cayley tree, e.g., for a spin chain (see 
below). Recall that there is no phase transition in spin 
models on a Cayley tree due to the presence of boundary 
spins. The number of these spins is of the order of the 
total number of spins on a given Cayley tree. 

In the ferromagnetic phase, i.e. at T < T c , or at H ^ 
the parameter x is smaller than 1. Note that in accor- 
dance with Eqs. I|21l) and lll'L'1) the non-local susceptibility 
x(i) of a regular Bethe lattice is finite at all temperatures 
and magnetic fields. 

Now let us consider an uncorrelated random network. 
We calculate a value of Xiji^ij — £) averaged over the 
ensemble of uncorrelated random graphs with a given 
degree distribution function P{k): 

X (i) = xa (4- = ^) = ££•••£ £ x« 



P(kj)kj J j-j- P{k n )k n {k n 1) | _P(fcj)fcj 



Zl 



,n=l 



"-=2 



(23) 



Here, the quantity in the square brackets is the proba- 
bility that a vertex i with degree fcj is connected with a 
vertex j having degree kj by a path that goes through 
vertices with degrees hi, &2, •••fef-i- 

In the paramagnetic phase at T > T c = 
2J/\n[(k 2 )/((k 2 } - 2{k))] and H = we have x tj = 1, 



and Eq. leads to the Eq. (|2*2*|) . Therefore, in zero 

magnetic field both on regular and random Bethe lat- 
tices the non-local susceptibility and Xij have the same 
temperature dependence determined by Eq. 122fl . 

At a magnetic field and in the ferromagnetic phase 
an approximate expression for the function xif) may be 
obtained in the framework of the approach proposed in 
|l2j |. We introduce Xij = exp(— hij). Here, the quantities 
hij are positive and independent random parameters. We 
use the following ansatz, 



E/kj » kh + Oik 1 ' 2 ), 



(24) 



where h = hij is the average value of the parameter. 
Applying the ansatz ffity to Eq. (fT5|l yields 

h = -i^P(k)k\ny{H,e-( k -^ h ) (25) 
k k 

which determines ft as a function of T and H. In fact, 
the parameter h plays the role of the order parameter. 
At H = 0, we get h = above T c and non-zero below T c . 
With this ansatz, 



where 



X(£)=z 2 1 z^- 1) AB^C, 



B = ^q(q-l)P{q) 
q 



( e H0 _|_ e -HP-qhy ' 

2 e -(8- 2 ) fe sinh(2J/3) 

( e (J+H)f3 + e -(J+H)/3-(q-l)h^2 

4 e -(q-i)h sinh (2 J/3) 
U7+H)0 + e -(j+H)0-( q -i)hy ■ 



(26) 
(27) 
(28) 
(29) 



The quantities A, B and C are finite (independent of N) 
for an arbitrary degree distribution P(k). The finitcness 
of the sums in Eqs. H27|) - (|29|l is ensured by the exponen- 
tial multiplier e~ qh . 

Result (|26f) is in agreement with the conclusions of the 
preceding section [compare relations l|26|l and ©]■ Note 
that the conclusion that x(£ > 2) cx G(£ > 2) — > as 
N — > ex) if the second moment of the degree distribution 
diverges, does not depend on the fact that we used the 
approximation i|24|). Indeed, substituting Eq. l(T%|l into 
Eq. I|23l) . one can prove that all the sums over degrees 
k n converge both in the ordered state and at H 7^ 0. 
Therefore, we get x(£) ^ z i ^ 1 ^ m agreement with Eq. 
(EH. 



IV. DISCUSSION AND CONCLUSIONS 

In the framework of a phenomenological approach and 
by using the Ising model we have studied time and static 
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pair correlations (xiffxjXjfo)) between interacting ob- 
jects clcfinccl on a random uncorrelatcd complex network. 
We have demonstrated that if a complex network has a 
divergent number of the second neighbors, then the pair 
correlations between second and more distant neighbors 
are strongly suppressed in large networks. In fact, the 
correlations tend to zero in the limit N — > oo. One should 
note that a particular form of the characteristic Z2(N) 
strongly depends on a network model. It is important 
that the size dependence of the pair correlation function 
is determined by the factor l/z^lN). 

The effect of weakening of the pair correlations is well 
known to people leaving in large towns and forming a 
large network of acquaintances. Close relationships are 
only being kept between close relatives inside family, 
close friends or between people working in the same of- 
fice. These people are our nearest neighbors. Our "sec- 
ond" and more distant nearest neighbors — people living 
on the next floor or working in a neighboring office - 
usually weakly influence our personal and social life. This 
is contrary to a small village where strong correlations are 
present between almost all inhabitants. 

Our analysis may be generalized to correlated net- 
works. Many natural networks exhibit correlations be- 
tween degrees of adjacent vertices, see, for example, Ref. 
pi). Note that in large uncorrelated networks, the diver- 
gence of the mean number zg of the £-th nearest neighbor 
can occur only simultaneously at all I > 2. In contrast, in 
large correlated networks, it is possible in principle that, 
say, 2:2,2:3 is finite and only z^>4 diverges. In this ex- 
ample, pair correlations are observable between the first, 
second and third nearest neighbors and vanish starting 



from the fourth nearest neighbors. 

We restricted ourselves to pair correlations. One 
should emphasize that pair correlations in cooperative 
systems on networks, unlike lattices, never show critical 
behavior. Even at the critical point of an interacting 
system on a network, we observe an exponentially rapid 
decay of pair correlations. In contrast, volume correla- 
tions in interacting systems on networks demonstrate a 
critical feature. Let us explain this point. The distri- 
bution of the full response of a system to a small lo- 
cal field is P(e) — J2i^(J2j Xij ~ e ); where 5(e) is the 
delta-function. This distribution is a rapidly decreasing 
function both below and above a phase transition on a 
network. On the other hand, at the critical point, P(e) 
is power-law, that is, "critical" [Toll22T|. 

One can conclude that a network structure of an in- 
teracting system strongly influences the dependence of 
pair correlations on the distance between interacting ob- 
jects. So, investigations of the correlations may be a use- 
ful method to understand the network topology. If an in- 
teracting system is defined on a network with a divergent 
number z n of n-th nearest neighbors of a vertex, then dy- 
namic and static pair correlation are strongly suppressed 
at distances n and greater. 
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